Equilibrium and dynamical properties of the ANNNI chain at the multiphase point 
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We study the equilibrium and dynamical properties of the ANNNI (axial next-nearest-neighbor Ising) chain at the multiphase 
point. An interesting property of the system is the macroscopic degeneracy of the ground state leading to finite zero-temperature 
entropy. In our equilibrium study we consider the effect of softening the spins. We show that the degeneracy of the ground 
state is lifted and there is a qualitative change in the low temperature behaviour of the system with a well defined low 
temperature peak of the specific heat that carries the thermodynamic "weight" of the ground state entropy. In our study 
of the dynamical properties, the stochastic Kawasaki dynamics is considered. The Fokker-Planck operator for the process 
corresponds to a quantum spin Hamiltonian similar to the Heisenberg ferromagnet but with constraints on allowed states. This 
leads to a number of differences in its properties which are obtained through exact numerical diagonalization, simulations and 
by obtaining various analytic bounds. 



PACS numbers: 05.50.+q, 02.50.Ey, 05.40.-a 



I. INTRODUCTION 

The ANNNI chain is one of the simplest systems with 
competing interactions. It is defined by the following 



Ising spin Hamiltonian: 



H 



i=l 



i+2) 



±1. 



(1) 



For J2 > 0, the interactions are competing and one can 
have different ground states depending on the relative 
strengths of the interactions. A specially interesting case 
is the point J\ — 2J 2 , the so-called multiphase point, 
where the ground state is no longer unique. It can be 
shown that any spin configuration, which does not have 
three consecutive spins of the same sign, is a ground state. 
For a chain of length L, the number of ground states 
~ fi L , where \i = (y/E + l)/2 is the golden mean. Thus 
there are an exponentially large number of degenerate 
ground states and the system has finite zero-temperature 
entropy per spin. The model has been extensively studied 
both in one and higher dimensions and is known to have 
a rich and interesting phase diagram Q . In this paper we 
consider some aspects of the equilibrium and dynamical 
behaviour of the ANNNI chain at the multiphase point. 

In our equilibrium study we consider the effect of soft- 
ening the spins, that is allowing them to take continuous 
instead of discrete values. It is usual in the study of 
spin models to consider soft-spin versions of discrete spin 
models. A well-known example is the Ginzberg-Landau 
Hamiltonian which is a continuum version of the discrete 
Ising model. Other examples occur in the study of spin 
glass models. For instance, the soft spin version of the 
Sherrington-Kirkpatrick (SK) [|| model was studied EJ 
in the context of dynamics. The reason for going to soft- 
spin versions is that they are often more amenable to 



theoretical approaches. It is usually expected that quali- 
tatively the soft and hard spin versions should show sim- 
ilar behaviour. 

For models with multiple ground states, arising out of 
frustration, softness may however change the degeneracy 
completely, as may be seen in a three-spin example, or 
as in the present case as we shall show here. The effect 
of spin-softening in systems with competing interactions 
has been studied earlier by several authors. Seno and 
Yeomans Q] have looked at the effect of softening spins 
at the multiphase point of a clock-model. They find, us- 
ing a perturbative method, that as a softness parameter 
is varied the system goes through a series of different 
ground states. In this work we use a similar perturbative 
method to prove that the macroscopic degeneracy of the 
ground state in the ANNNI model is lifted by the small- 
est amount of softness. We then show explicitly how the 
release of the zero-temperature entropy results in qual- 
itative differences in the low temperature properties of 
the system. This is similar to the recently observed phe- 
nomena of entropy release in spin- ice systems ||. We 
also construct an effective hard-spin Hamiltonian to de- 
scribe the low- temperature properties of the soft-spin 
model. We have also performed Monte Carlo simulations 
on the soft-spin model and verified the low-temperature 
predictions of the effective Hamiltonian. 

In the second part of the paper we look at the dy- 
namical properties of the system. As noted before, the 
ANNNI model at the multiphase point has a large num- 
ber of degenerate ground states. It is, therefore, of in- 
terest to look at dynamical properties of the system at 
low temperatures. Here we use Kawasaki dynamics to 
evolve the system and consider zero temperature prop- 
erties only. Thus two nearest neighbor spins flip with 
a rate 7, provided both magnetization and energy are 
conserved. This dynamics has been studied earlier by 
Das and Barma In this paper we extend their stud- 
ies by using the correspondence between VF-matrices for 
stochastic processes and quantum spin chains. 

The correspondence between the stochastic Fokker- 
Planck operator and quantum chains has often been ex- 
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ploited to derive dynamical properties. For instance, the 
scaling, with system size, of the first excited state of the 
quantum Hamiltonian gives the dynamical exponent of 
the stochastic process. A well-known example where this 
correspondence has been used is in exclusion processes 
0, which are stochastic models of hard-core diffusing 
particles. For such processes, it has been possible to 
exactly calculate the dynamic exponent by solving the 
corresponding quantum model, namely the Heisenberg 
model ||. The dynamics considered by us is very sim- 
ilar to the symmetric exclusion process (SEP) but with 
added restrictions on allowed moves. To see this we first 
note that with the present dynamics, nearest neighbor 
spins with opposite signs flip, provided that the resulting 
configuration satisfies the ground state constraint of no 
three successive spins having same signs. Identifying up 
spins with particles and down spins with holes we see that 
the dynamics is equivalent to hard core particles diffus- 
ing on a lattice with the constraint that there cannot be 
three successive particles or holes. An interesting ques- 
tion is whether these rather strong constraints make the 
system different from the SEP. Earlier numerical work 
H seems to suggest that the dynamics still behaves like 
the SEP. We note that there have been some other re- 
cent studies on exclusion processes with constraints on 
allowed configurations These cases are solvable by 
the Bethe ansatz and show the same behaviour as the 
unconstrained model. 

Here we address this question of the effect of the 
constraints by studying the quantum Hamiltonian. By 
means of exact numerical diagonalization for finite chains 
and through analytic bounds, we have tried to un- 
derstand the differences and similarities of the present 
Hamiltonian with the Heisenberg Hamiltonian for the 
SEP. We also discuss the different symmetry properties 
of the two quantum models. The Heisenberg model has 
full rotational symmetry and this has several important 
implications some of which are of direct relevance in un- 
derstanding the original stochastic process. For example 
it implies that two-point time correlations in the SEP 
do not depend on the number of particles. The present 
model, on the other hand is only invariant under rota- 
tions in the XY plane. 

The rest of the paper is divided into two sections. In 
section (II), we consider equilibrium properties of the 
soft-spin model while in section (III) we consider the dy- 
namics of the hard-spin model. Section (IV) contains a 
summary of our main results and a few concluding re- 
marks. 



II. SOFT-SPIN ANNNI MODEL 

We consider the following soft spin version of the 
ANNNI model: 

H S = Y^ J{2s iSl+1 + SiS i+2 ) + a g(sf/i - sf/2), 



Si S (-00, 00) (2) 

where a is a dimensionless parameter which controls the 
amount of softness. In the limit a — > 00 we get the hard- 
spin model. We will set g = 1 since there is no loss of 
generality in doing so. 

Let us first look at the ground states of the soft-spin 
Hamiltonian given by Eqn. (|J). To do so we look at the 
extrema of H s which are obtained by setting dH/dsi = 
for all i. This gives: 

2J(s i+ i + Si-i) + J{s l+2 + Si_ 2 ) + a(sf - Si) — (3) 

Solving this set of coupled nonlinear equations in general 
is very difficult. However for small values of the param- 
eter I/a we can obtain the solutions perturbatively. For 
a — ► 00 all configurations, {s^}, with s$ = 0,±1 are 
solutions. Those with {si = ±1} correspond to the min- 
ima. For finite but large a we try to obtain the solutions 
perturbatively with 1/a acting as the perturbation pa- 
rameter. We denote the unperturbed minima by the set 
{U = ±1}. Let us try the following perturbative expan- 
sion: 

00 -. 

where the coefficients ij are independent of a and = 
ti = ±1 correspond to the unperturbed solutions in the 
limit a — ► 00. Plugging this into Eqn. ([?]), we get 

J[2{t l+l + U-i) + (U+2 + ti_a)] + -[2(4+ } i + t£.\) + 

a 

(*& + *&)] + *P + |(3t,(*^) a + 2t?>) + 0(1) = 0. 
Equating different powers of 1/a to zero we then get: 

if 5 = ^-[2(t i+ i+ti_ 1 ) + (t i+2 +t i _ 2 )]+0(i) 

* (2) =^[2(e+e\)+(e+e ) 2 )]-^(f ) ) 2 

and so on. Thus we get 2 L perturbed minima given by 
the above perturbation series. The energies correspond- 
ing to these minima can now be found by putting these 
solutions into the expression for energy in Eqn. (0). We 
thus get 

E = E + Ex+ E 2 + 0{\/a 2 ), where 




E\ — J (2ij£i+i + titi+2) 

i 

E * = ~ Et(^ (1) ) 2 + 2J M+1 + *i+l*f 3 ) + J (Ut\% + ti+2-t^)} 

a * — ' 

i 

_J2 

= -7T- E( 5 + 4< ^+l + 4 *i*i+2 + 4^+3 + tiU+i). (5) 
la z — ' 
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In the above expansion, Eq corresponds to the unper- 
turbed energy, while E\ and E2 represent the corrections 
resulting from the perturbation. In the a —* 00 limit the 
term E\ causes the energy levels of the 2 L minima to 
split, with separation between them ~ O(J). We rec- 
ognize Ei as the Hamiltonian for the hard-spin ANNNI 
model. Thus the lowest energy level is still /i L -fold de- 
generate. The term E2 then causes a further splitting of 
the ground states into levels with separation ~ 0{J 2 /a). 

To see whether or not the macroscopic degeneracy of 
the ground state survives, we need to consider the inter- 
action Hamiltonian corresponding to the energy term Ei . 
Since we are interested in the splitting of the lowest en- 
ergy level of E\ , we only consider the restricted subspace 
of spin configurations which are ground states of E\. In 
this subspace the Hamiltonian corresponding to Ei can 
be rewritten as 



-3LJ 2 J 2 



2 a 



2a 



+ 4^+3 + Uti+i). (6) 



Thus all the interactions are ferromagnetic. However the 
ground state of H2 is not the state with all spins up, since 
this does not belong to the subspace of ground states of 
Ei. To find the ground state, we write the second term 
in H2, which we denote by /12, in the following form (the 
constant factor J 2 /(2a) is suppressed): 

hi = — ^^(2tit i+ 2 + 4£j£j+3 + titi + 4) 

i 

= — J]] e(ii,*i+l)ij+2,*i+3 I ti+4,iti+5,ti+6)ti+r) 
i=(4«.+l) 

where e(ti, t 2 , t 3 , t 4 | t 5 , t 6 , t 7 , t 8 ) = tit 3 + t 2 t 4 + 2t 3 t 5 + 
2* 4 * 6 + t 5 t 7 + * 6 * 8 + 2*1*4 + 4*2^5 + 4t 3 * 6 + 4t 4 ^ + 2t 5 t$ 
+*i* s + t 2 t 6 + * 3 *7 + * 4 *8 (7) 

and the index n runs from to (L/4 — 1) (we take L to 
be an integral multiple of 4). By enumerating the matrix 
elements e(*i, *2, £3, *4 | t5,te,t 7l ts) for all allowed spin 
configurations we find that the lowest energy configura- 
tion is obtained for the periodic sequence (Tt-lttl •■■) an d 
the five other configurations obtained by translating and 
flipping this. Thus we find that the infinite degeneracy of 
the ground state is removed and instead we get a six-fold 
degenerate ground state. We note that the procedure just 
outlined provides a straight forward method of finding 
the ground state of any spin Hamiltonian. By numeri- 
cally solving Eqn. (S) for small lattice sizes (L — 12) and 
finding the minimum energy configurations for a large 
enough (a = 50), we have verified that the perturbative 
solutions are quite accurate. 

The fact that softening of the spins results in removal 
of the exponential degeneracy of the ground state means 
that the finite zero temperature entropy is released and 
we expect it to show up in the behaviour of the low 
temperature specific heat. This leads to the soft-spin 
model having low-temperature properties very different 
from the hard-spin version as we shall now see. 



We note that the hard spin model is easily solvable by 
transfer-matrix methods and one can exactly compute 
various thermodynamic properties. In the soft-spin case 
the transfer-matrix eigenvalue equation becomes an in- 
tegral equation which we have not been able to solve. 
Hence we have studied the model by Monte Carlo sim- 
ulations. We have used a dynamics which allows three 
kinds of processes; 

(i) single spin-flip moves, 

(ii) moves in which two nearest neighbor spins are si- 
multaneously flipped and 

(iii) moves which change the length of a spin. 

All three kinds of processes occur with usual Metropo- 
lis rates. The reason for allowing both single and double 
spin-flips is the following. We find that in the hard-spin 
case, equilibration times, with a single-spin flip dynam- 
ics, become very large at low temperatures. On the other 
hand, allowing for two-spin flips results in very fast equi- 
libration. This is related to the fact that while the single 
spin-flip dynamics at T = is non-ergodic, including 
double- flips makes it ergodic. We expect a similar situ- 
ation even in the case of soft-spins and so have included 
both (i) and (ii) . Finally (iii) is necessary since the spins 
are now continuous variables and we need to be able to 
change their lengths. 

In order to compare the properties of the soft-spin 
model with those of the hard-spin one, it is necessary 
to subtract from the soft-spin free energy a part corre- 
sponding to the continuum degrees of freedom. We thus 
look at the following free energy: 



F = (-l//3)[lnTre- 0Hs +iln(2) - InTre 



(8) 



where H s is as in Eqn. H g = £\ a ( s t /4 — sf/2), and 
Tr indicates integration over all spin variables. We note 
that the above expression for the free energy is equivalent 
to writing the partition function in the form 



p ( s_ )=n 



Z = Tre- fJH P(s) with 

2e -/3a( S ?/4- s =/2) 



< J ds * 



-0a(«f/4-«J/2) ' 



(9) 



H being the original hard-spin Hamiltonian and P(s) a 
probability distribution over the spin variables. In the 
limit a — > 00 this exactly reduces to the hard-spin parti- 
tion function while at T — > 00 one gets Z = 2 L . ^From 
our simulations we get properties corresponding to the 
first part of the free energy in Eqn. The second 

part simply corresponds to a nonintcracting system and 
its properties can be easily computed numerically. 

In Fig. [j]we plot the specific heat data C(T) for both 
the soft-spin and hard-spin models. The hard-spin result 
is exact and corresponds to infinite system size while the 
soft-spin data is from simulations on a chain of length 
L = 24. The values of various parameters used in the 
simulation were a — 50 and J = 1. The high temperature 
(T > 1) data was obtained by averaging over 10 6 Monte 
Carlo steps while the low temperature data is over 10 7 
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steps. As expected we find a second peak in the specific 
heat at low temperatures. For the hard-spin case the 
total area under the curve for C(T)/T is equal to ln(2//i). 
The ground state entropy, ln(/i), which is released when 
the spins are softened, is mostly accounted for by the 
area under the low temperature peak. 

The low temperature properties are quite well repro- 
duced by the effective Hamiltonian, H2, which describes 
the energy levels in the lowest band. The thermodynamic 
properties of H2 can be exactly calculated by transfer 
matrix methods, both for finite system sizes and in the 
infinite size limit. In Fig. ^ we plot the soft-spin low 
temperature simulation data C(T) for two system sizes 
and compare them with results obtained from the effec- 
tive Hamiltonian. We see good agreement between the 
two. We also show the infinite system size C(T) curve 
obtained from the Hamiltonian Hi- It is interesting to 
note that the peak value of the specific heat first increases 
with system size and then starts decreasing beyond a cer- 
tain size. 



III. KAWASAKI DYNAMICS OF THE 
HARD-SPIN ANNNI MODEL AT THE 
MULTIPHASE POINT 

As for the usual exclusion process, the quantum Hamil- 
tonian corresponding to our process can be easily written 
and is given by: 



1 



k=l 



y k °i+i-i)]p k }v (10) 

where crj? are the usual Pauli matrices, Pk are local pro- 
jection operators given by 



P fe = (1 - 07c-207c-i)(l - erfc +2 cr fc+ 3)/4 



(11) 



and V = H fe=1 Pk is a global projection operator which 
projects onto the space of allowed states, i.e those that 
satisfy the ground state constraint. The spin- flip rate, 
7, has been set to unity. Alternatively we can write the 
Fokker-Planck operator in the following form: 



n = -J2( e k + °l) where 



k=l 



(12) 



The term J^k ®k ^ s l ne diagonal term since it corresponds 
to flipping an unequal pair twice. It is important to write 
the diagonal part carefully. For instance if in Eqn. (10), 
the local projection operators, Pk, were not present, the 
off-diagonal elements of Tt would still be correct but the 
diagonal ones would be wrong. 



We now study the properties of this quantum Hamil- 
tonian. Our interests are (a) to compare the symmetry 
properties and conservation laws of the present Hamilto- 
nian with that of the Heisenberg model and (b) to obtain 
results on the energy gap and hence the dynamical expo- 
nent. 



A. Symmetry properties and conservation laws of 
the quantum model 

We first observe that the z-component of the total spin, 
S z , commutes with TL. This simply implies conservation 
of spin or number of particles in the stochastic model. 
Thus we can classify energy states into sectors labelled by 
number of particles n. The constraints on allowed config- 
urations means that for a lattice of length L the number 
of particles can vary over the range [L/3] < n < L— [L/3] 
where [L/3] denotes the smallest integer greater than or 
equal to L/3. It can be shown that except in the lowest 
and highest sectors, in every other case the dynamics is 
ergodic. It then follows from detailed balance that the 
steady state is one in which all allowed configurations 
in a given sector occur with equal probability. For the 
quantum model this means that the ground state in any 
sector is an equally weighted sum over all states ( For 
the special case where L is a multiple of 3, the lowest 
and highest sectors have 3-fold degenerate ground states 

The other components of the total angular momen- 
tum S x and S v however do not commute with Ti.. Thus 
the present Hamiltonian has U(l) symmetry instead of 
the SU(2) symmetry of the Heisenberg model. Also 
even though the ground states are degenerate, with one 
state in every S z sector, there is no analogue of the rais- 
ing/lowering operator . If there were such an operator 
then the entire eigenvalue spectrum in the n-particle sec- 
tor would be a subset of the (n — l)-particle sector (for 
n < L/2). By looking at the spectrum for finite sized 
lattices we have verified that this is not so. 

To study the presence of long-range order in the 
ground-state, we have calculated the two-point static 
correlation functions c z (r) =< a^a z > and c±(r) =< 
0o~<7~ > in the ground state for the half-filled sector. The 
simple characterization of the ground states in terms of 
disallowed subsequences enables calculation of ground- 
state expectation of any operator by means of trans- 
fer matrices. The transfer matrix method sums over 
all the different particle sectors, but in the thermody- 
namic limit the half filled sector dominates, and so we get 
correct results (To compute expectation values in other 
sectors one would need to introduce a chemical poten- 
tial). Thus we find that c z (r) = Acos(cj) - 2irr/3) e - r ^ 
where £ = l/log((3 + v / (5))/2) = 1.03904... and A and 
<p are constants that have different values on odd and 
even sites. Fourier transforming c(r) gives the structure 
factor < a z \—q)a z (q) > which has the form shown in 
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Fig. H We note that it is non-vanishing at all q. The 
off-diagonal correlation can similarly be obtained using 
transfer matrices but the calculation becomes extremely 
cumbersome. Instead we have computed this correla- 
tion numerically for finite lattices and find that it satu- 
rates, for large r, to a constant value, which is given by 
< <t_ > 2 = 0.02917... (which has been obtained by using 
the transfer matrix method). 

Thus we find that ground-state correlation functions 
show the same behaviour as in the Heisenberg chain. 
For the Heisenberg model c z (r) is delta-correlated while 
c±(r) saturates to the value 1/4 (which is much larger 
than its value in the present model). The presence of 
off-diagonal long range order means that U(l) symme- 
try is broken in the ground state. This is analogous to 
the breaking of SU (2) symmetry in the ground state of 
the Heisenberg model. On the other hand, consider the 
XX Z chain |l(J defined by the Hamiltonian 

L A 

k=l 

Away from the two isotropic points (A = ±1), this has 
the same symmetry as the present model. It has no long 
range order in the gapless phase (— 1 < A < 1) and all 
correlations < a a (0)cr a (r) > have power law decays. In 
the ferromagnetic phase (A > 1), the model has a gap 
and full ferrromagnetic long-range order in the ground 
state, with ultra-local longitudinal correlations namely 
c 2 (r) = 1/4. Thus we see that as far as ground state 
correlations are concerned the present model is different 
from the anisotropic XX Z chain even though they have 
the same symmetry properties. Our model is more sim- 
ilar in properties to the ferromagnet (A = 1) but has a 
nontrivial depletion of the condensate, as well as a non- 
trivial < a^a z > correlation. 

Finally we note that rotational invariance of the 
Heisenberg model means that two-point time correlations 
are completely determined by single magnon excitations 
and so have the same behaviour in any S z sector |.|n). 
This result does not hold in the case of the present model. 

A second conserved quantity in the model is the total 
linear momentum. This follows from the translation in- 
variance of TL. The momentum operator commutes both 
with TL and S z so that in each S z sector energy states 
can be labelled by their momentum. Clearly the ground 
state has zero momentum. 



B. Results on the energy gap 

As is well known the first excited state of TL deter- 
mines the decay of correlations for the stochastic process. 
Thus the energy gap A ~ 1/L Z and this determines the 
dynamic exponent z. For the SEP, which corresponds 
to the Heisenberg ferromagnet, it is known that z = 2. 
This simply reflects the diffusive modes in the dynamics. 



The dynamics studied here is very similar to the SEP but 
with the constraints on the allowed number of succesive 
particles and holes. An interesting question is whether 
these rather strong constraints change the dynamical ex- 
ponent. Unlike the Heisenberg model where the Bethe 
ansatz is applicable and yields information on the eigen- 
value spectrum, the Hamiltonian in Eqn. ( )l0| ) is much 
more complicated and we have not been able to use the 
Bethe ansatz. We have looked at the eigenvalue spec- 
trum by numerical diagonalization of TL for small system 
sizes and also through Monte Carlo simulations. We also 
obtain various analytic bounds on the energy levels. 

(i) Results of numerical diagonalization of TL and 
Monte Carlo simulations 

We have carried out exact diagonalization of the 
Hamiltonian in Eqn. ( |ici| ) for chains of length upto 
L = 22 at half filling. The diagonalization has been done 
in the momentum basis. This makes the Hamiltonian 
block diagonal and enables us to go to quite large chian 
sizes. We find that for small L the first excited state oc- 
curs at total linear momentum q — ix and the gap seems 
to decreases as ~ 1/L. However from L = 22 onwards, 
the first excited state shifts to q = 2n/L and the gap 
at this momentum decreases as ~ 1/L 2 . In Fig. [| we 
show the numerically obtained gaps at the two momenta 
as a function of system size. We also plot correspond- 
ing upper bounds on the gaps (to be derived in the next 
section). 

We note here that though it is usually the first ex- 
cited state that determines the decay of correlations in 
the stochastic process, it is possible to construct corre- 
lation functions whose decay is governed by some other 
eigenvalue. As an example consider the operator Q — 
e l ~12k k ' Jk . This is the so-called twist operator, first 
studied by Lieb, Schultz and Mattis ]l2] |. In this case, the 
decay of the correlation, < Q(0)Q(t) > is determined by 
the lowest eigenvalue at momentum it since the operator 
carries momentum tt. In Fig. |^ we show the decay con- 
stant as determined from the correlation decay for differ- 
ent system sizes and compare them with those obtained 
from exact diagonalization. The correlation function is 
obtained from Monte Carlo simulations and can also be 
used for larger system sizes at which numerical diagonal- 
ization becomes too difficult. 

(ii) Exact Bounds 

We now find upper bounds on the first excited state. 
Consider the sector with states which have n overturned 
spins. The bounds are obtained by constructing trial 
wave functions orthogonal to the ground state in each sec- 
tor. Thus consider the operators o- z (q) = -^^Zk a k elkq 
and the twist operator Q defined in the previous section. 
Under translation these operators transform as 
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Ta*(g)Tt=-L£a fe+1 , 
TQT^ = e^S fc fcCT J+i 



ifcg 



V( g ) 



e i2,rd Q 



where <i = u/£ is the filling fraction of particles. If 
|0„ > is the ground state in the n-particle sector, then the 
states <7 2 (<7)|0„ > and Q|0„ > have momenta q and 27rd 
respectively and for q ^ are orthogonal to the ground 
state, which has zero momentum. Hence the following 
expectation values give us two different upper bounds on 
the gap: 



(a) 
(b) 



_ < <T z (-q)H<r z (q) > 
< o- z {-q)a z {q) > 
= < Q^HQ > 



(14) 
(15) 



where < ... > denotes ground state expectations. We 
now evaluate (a) and (b). We shall henceforth restrict 
ourselves to the half-filled sector only, though extensions 
to other sectors can be done. 

(a) To evaluate e z we first note that the numerator and 
denominator in Eqn. ( jl4| ) can be written in the following 
equivalent form: 

< a z (-q)Ha z (q) > = \ £ < [a z , [H,vf +1 ]] > 

i 

< a z (-q)a* e< * < ^°f+l > ( 16 ) 



The commutator occurring in the above equation can be 
evaluated and gives: 

[al [H,of +1 ]] = -4r[(a+az + ^a+) Pl (-d hL + Si A ) + 
4(<7+^- + o- L ot)vL{^L-i - Kl)\V^ (17) 



Inserting this in Eqn. (|16| ), and using translational in- 
variance of the ground state we finally obtain: 

<a z (-q)Ha z (q) > 

= 4[1 - cos(g)] < VPi(<t+<t2 + Ox a t)^ > 
= 2[l-oos(g)] <VP 1 {\-a\a z 2 )V> (18) 

where the last step has been obtained using the fact that 
< 0|H|0 >= 0. As noted before ground-state expec- 
tations of any operator can be computed using transfer 
matrices. The expectation value on the rhs of Eqn. jlq ) 
is thus found to have the limiting value (as L — ► oo)< 
VPi(l-a(a$)V >= 8-16/ VS. The Fourier transform of 
c z (r), which gives the structure factor < a z (—q)a z (q) >, 
has already been obtained and was plotted in Fig. 0. We 
note that it is non- vanishing at all q. Finally, from Eqns. 
( p^ , ^6| ) we get e z which is plotted in Fig. || along with 
the exact results from finite size diagonalization. Putting 
q = 2n/L and putting in all numerical factors, we get the 
following result: 



(b) We now obtain the other bound using the twist 
operator, Q. We first note the following properties of Q: 

Q^a+<j- +1 Q\{<i} > = e l ^a+a w \{a} > 
Q^afa+ +1 Q\{a} > = a+ +1 \{a} > . (20) 

Using these relations we obtain 

< 0\Q^HQ\0 > 

k 

- E \ < o|O^K^. +1 - l)P k VQ\0 > 

k 

= -cos(27r/L)]r < 0\V(a+<j- +1 +a-<j+ +1 )P k V\0 > 

k 

-Y,\<v\n°ivi + x-VPkV\u> 

k 

= |[1 - cos(27r/i)] < Q\V(1 - cj z k al +1 )P k V\0 >, (21) 

where in the last step we have again used < O|7i|0 >= 
and translational invariance of the ground state. The ex- 
pectation value above has already obtained so that we 
get, for large L, the following bound for the gap at mo- 
mentum q = tt. 



A < 0.845— 
L 



(22) 



A < 19.78— r 



(19) 



In Fig. |J we have plotted both the bounds and the exact 
finite size results at q = 2ir/L and q = tt as functions of 
the system size. 



IV. SUMMARY 

In summary we have studied a one-dimensional spin 
model with competing interactions and studied its low- 
temperature equilibrium and dynamical properties. In 
the equilibrium case we have shown that low tempera- 
ture properties of the soft-spin and hard-spin versions of 
the model can be very different. The hard-spin version 
of the model has an infinitely degenerate ground-state. 
Through a perturbative calculation we have shown, that 
as soon as we introduce the slightest amount of softness, 
the degeneracy is lifted. The ground state energy levels 
split to form a band which is separated from higher lev- 
els by AE — O(J). The energy levels within this lowest 
band are described by an effective hard-spin Hamilto- 
nian, containing ferromagnetic interactions upto fourth 
neighbour terms. This can be used to approximately de- 
rive the low temperature properties of the model. Wc 
find reasonably good agreement with results from Monte 
Carlo simulations of the soft-spin model. 

Our results indicate that the fixed-length (a — > oo) 
limit is a singular one in our model at low temperatures. 
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Since the ground state of the soft-spin model for large but 
finite a is only six-fold degenerate, it would order into one 
of these six states as T — > 0. This implies the occurrence 
of a zero-temperature phase transition and the existence 
of an appropriately defined correlation length that di- 
verges as T goes to zero. In the fixed-length (a — > oo) 
limit, on the other hand, averaging over all the degener- 
ate ground states leads to a finite correlation length even 
at T = 0. These results suggest that it would be inter- 
esting to study the effects of softening the spins on the 
thermodynamic behavior of two- and higher-dimensional 
hard-spin models with extensive ground-state entropy. A 
well-known model of this kind is the nearest-neighbor 
Ising antiferromagnet on a triangular lattice This 
model does not exhibit any phase transition at a non-zero 
temperature. The degeneracy-lifting effect of introducing 
magnitude fluctuations found in our study suggests that 
soft-spin versions of this and other similar models may 
exhibit finite-temperature phase transitions. Further in- 
vestigation of this question would be very interesting. 

We believe that the removal of the exponential ground- 
state degeneracy by the introduction of spin-softness in 
the model studied here is a special case of a more gen- 
eral phenomenon in which the presence of additional de- 
grees of freedom allows the system to relieve frustra- 
tion and thus reduce the number of degenerate ground 
states. Coupling the hard spins to other degrees of free- 
dom, such as elastic variables describing possible defor- 
mations of the underlying lattice, would probably have 
similar effects on the degeneracy of the ground state. It 
is interesting to note in this context that a "deformable" 
Ising antiferromagnet on a triangular lattice in which the 
Ising spins are coupled to elastic degrees of freedom ex- 
hibits a Peierls-type phase transition at a non-zero 
temperature. The ordering of the spins at this transition 
is accompanied by a distortion of the lattice. In general, 
it is expected that in real, physical systems, such cou- 
plings to other degrees of freedom, however weak, would 
induce some kind of ordering of the spins as the tem- 
perature is reduced toward zero, thereby avoiding the 
unstable situation of having a non- vanishing entropy per 
spin at T = 0. 

Many disordered spin systems, such as spin glasses [[l6) , 
exhibit a large number of nearly-degenerate metastable 
states arising out of frustration. To take an example, 
the SK model of infinite-range Ising spin glass is 
known [l6| to have an exponentially large number of lo- 
cal minima of the free energy (locally stable solutions 
of the TAP equations (l5)) at sufficiently low temper- 
atures. These local minima of the free energy become 
local minima of the energy at T = 0. The presence of a 
large number (divergent in the thermodynamic limit) of 
nearly-degenerate metastable states is crucial in the de- 
velopment of the present understanding fl6| of the equi- 
librium and dynamic properties of this system at low 
temperatures. Our results about the lifting of degen- 
eracy by the introduction of spin-softness raise the fol- 
lowing interesting question: would the low-temperature 



properties of a soft-spin version of the SK model differ 
in any significant way from those of the original model? 
While soft -spin versions of the SK model have been used 
in studies p] of the dynamics, questions about how the 
number and properties of the metastable states of this 
model change as the spins are made soft have not been 
addressed in detail. Further investigation of these issues 
would be most interesting. 

Finally it is interesting to note that a similar way of 
lowering frustration is to make the coupling constants 
soft while keeping the spins hard. For example, in the 
case of the Edwards- Anderson Ising spin glass model, two 
versions have been studied fP7| ]. One is the ±J model 
where the nearest-neighbor coupling constants randomly 
take the discrete values ±J with equal probability. In 
the other case, the Js are chosen from a gaussian distri- 
bution. In d = 2, both these cases are believed to have 
zero-temperature phase transitions, but the nature of the 
transition is different in the two cases. This difference 
again arises because of the different ground-state dege- 
naracies in the two cases. In the ±J model, the ground- 
state is exponentially degenerate, while it is unique (mod- 
ulo a global inversion of all the spins) in the gaussian case. 
However in higher dimensions where the transition tem- 
perature is finite, critical properties near the transition 
appear to be the same in both cases. 

In our nonequilibrium studies we considered the 
Kawasaki dynamics and studied the quantum Hamilto- 
nian corresponding to the Fokker-Planck operator for the 
stochastic process. The spectrum of the Hamiltonian is 
obtained by numerical diagonalization of finite chains. 
An interesting crossover of the first excited state from 
momentum n to 2n/L is observed with increase in system 
size. We have found analytic upper bounds on the gaps 
at these two momenta. These, along with our numeri- 
cal diagonalization results suggest that the gap vanishes 
as ~ 1/L 2 and so the dynamics is diffusive as in SEP. 
We have also compared the symmetry properties of our 
Hamiltonian with the Heisenberg model. We find that 
while the model has the symmetry of the XXZ model, 
its ground-state properties are closer to those of the fer- 
romagnetic isotropic point. In summary we have shown 
that our model is a very nontrivial cousin of the Heisen- 
berg ferromagnet. The exclusion of three adjacent like 
spins essentially changes the model dynamics, and re- 
sults in a nontrivially depleted condensate in < cr^cT* > 
and a nontrivial gapped < u^al > correlation function. 
The existence of a groundstate in every S z sector is quite 
obvious from the stochastic point of view but a nontriv- 
ial one within the framework of the quantum system (e.g. 
the absence of a a~ operator), and require a deeper un- 
derstanding. 



7 



V. ACKNOWLEDGEMENTS 



We thank Mustansir Barma, Dibyendu Das and 
Deepak Dhar for useful discussions. AD is grateful to 
the Poornaprajna Institute for financial support. 



J. M. Yeomans, Solid State Physics 41, 151 (1988). 

D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 

1972 (1975). 

H. Sompolinsky and A. Zippelius, Phys. Rev B 25, 6860 
(1982). 

F. Seno and J. M. Yeomans, Phys. Rev. B 50, 10385 
(1994). 

A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, 

B. S. Shastry, Nature 399, 333 (1999). 

D. Das and M. Barma, Phys. Rev. E 60, 2577 (1999). 
T. Ligget, Interacting Particle Systems, New York, 
Springer (1985). 

L. H. Gwa and H. Spohn, Phys. Rev. A 46, 844 (1992). 
F. C. Alcaraz and R. Z. Bariev, Phys. Rev. E 60, 79 
(1999). 

C. N. Yang and C. P. Yang, Phys. Rev. 147, 303 ( 1966) . 
M. D. Gryn berg and R. B. Stinchcombe, cond 



mat/9908051/ 



E. Lieb, T. Schultz and D. Mattis, Annals of Phys. 16, 
407 (1961). 

G. H. Wannier, Phys Rev. 79, 357 (1950); R. M. F. 
Houtappel, Physica 16, 425 (1950). 

Z. Y. Chen and M. Kardar, J. Phys. C 19, 6825 (1986); 
Lei Gu el al, Phys. Rev. B 53, 11985 (1996). 
D. J. Thouless, P. W. Anderson and R. G. Palmer, Phil. 
Mag. 35, 593 (1977). 

M. Mezard, G. Parisi and M. A. Virasoro, Spin glass 
theory and beyond (World Scientific, Singapore, 1987). 
K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 
(1986). 



- soft spins 
hard spin 




Fig. 1 Dhar, Shastry, Dasgupta 



FIG. 1. Simulation data C(T) for the soft-spin model on a 
lattice of size N = 24. A low temperature peak can be seen. 
For comparision we have also plotted the hard-spin results. 
Most of the entropy released (~ 85%) is contained within the 
low temperature peak while the remaining occurs in the high 
temperature region (shaded portion). 



Monte Carlo: N=12 
Monte Carlo: N=24 
Perturbative: N=12, 24 
Perturbative: N=oo 




Fig. 2 Dhar, Shastry, Dasgupta 



FIG. 2. The plot of C(T) at low temperatures as obtained 
from simulations and from the effective Hamiltonian for dif- 
ferent system sizes. We also show the effective Hamiltonian 
result for infinite system size. 
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Fig. 3 Dhar, Shastry, Dasgupta 



0.02 



0.04 



0.06 

1/L 



0.08 



0.1 



Fig. 5 Dhar, Shastry, Dasgupta 



FIG. 3. The diagonal structure factor plotted as a function 
of the total wave number q. 
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Fig. 4 Dhar, Shastry, Dasgupta 

FIG. 4. In this figure the exact energy gaps A at the two 
momenta 7r and 2n /L are plotted against inverse system size. 
Also plotted are exact bounds at the two momenta. 



FIG. 5. The energy gap, as obtained from the decay of 
the correlation function < Q(0)Q(t) > is plotted as a func- 
tion of inverse system size. Also plotted are the results from 
exact numerical diagonalization and the upper bound. The 
diagonalization has been done till system size L = 22 while 
the < Q(0)Q(t) > data is from Monte Carlo simulations for 
system size upto L = 36. 




Fig. 6 Dhar, Shastry, Dasgupta 



FIG. 6. The gap upper bound, e z , plotted as a function 
of total momentum, q. The exact eigenvalues for a system of 
size L — 18 are also shown. 
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